February 2, 2008 9:11 WSPC/Trim Size: 9.75in x 6.5in for Proceedings 



Main 



STABILITY OF ACCRETION DISK AROUND ROTATING BLACK 

HOLES 



BANIBRATA MUKHOPADHYAY * 
Astronomy Division, P.O.Box 3000, University of Oulu, FIN-90014, Finland 



I discuss the stability of accretion disks when the black hole is considered to be rotating. 
I show, how the fluid properties get changed for different choices of angular momentum of 
black holes. I treat the problem in pseudo-Newtonian approach with a suitable potential 
from Kerr geometry. When the angular momentum of a black hole is considered to 
be significant, the valid disk parameter region affects and a disk may become unstable. 
Also the possibility of shock in an accretion disk around rotating black holes is checked. 
When the black hole is chosen to be rotating, the sonic locations of the accretion disk 
get shifted or disappear, making the disk unstable by means of loosing entropy. To bring 
the disk in a stable situation, the angular momentum of the accreting matter has to be 
reduced/enhanced (for co/counter-rotating disk) by means of some physical process. 



The accretion phenomena around black holes is indeed a topic extensively dis- 
cussed over the last three decades (e.g. Mukhopadhyajd and references therein). 
Here I like to concentrate on the effect of black hole's rotation to the fundamental 
disk properties. I particularly show how the fluid properties get changed with the 
variation of black hole's angular momentum. Therefore the basic characteristics of 
a disk along with its stability around Schwarzschild black hole may differ from that 
of a Kerr black hole, always quantitatively and in some cases qualitatively also. 

Here the word stability has a very specific meaning. We know from the second 
law of thermodynamics, entropy is the measure of disorderedness of a system. As 
the disorderedness rises up, entropy increases, the system becomes more stable, and 
vice versa. With the change of fluid parameters, when disks gain or loose entropy 
at a certain region then the stability of that region is thought to be increased or 
decreased respectively. 

I will pursue the present work in the pseudo-Newtonian approach bringing a 
pseudo-potential applicable for a disk around Kerr black hole. More particularly, 
I like to use an effective gravitational force corresponding to a pseudo-Newtonian 
potential. Following an earlier worlpl such a gravitational force was proposed^ 
by me, which will be used here. An interesting point to note that this was a 
potential established for first time directly from the space-time metric. In that 
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respect I proposed an algorithm also, which has already been used in other works 
(e.g. GhosrP). 

Here the force corresponding to pseudo-potential is given as 

y ' x 3 (y/x(x - 2) + a) 2 w 

where a is specific angular momentum of a black hole and x is disk radius in unit 
of GM/c 2 , M is mass of the black hole. Throughout this paper, I also express the 
specific angular momentum and velocity in unit of GM/c and c respectively; others 
have their usual meaning. This F(x) will be used in radial momentum balance 
equation to incorporate relativistic effects approximately. 

As my main interest is to check the sole effect of black hole's rotation on the 
disk properties, here I consider the inviscid fluid so that there is no dissipative force 
into the system and the angular momentum of disk fluid is constant for a particular 
case. The set of model equations for a disk is given as 

d , , s s dv 1 dP A 2 

Tx t*ph( X )v)=0; v- + ---^ + F( X )=0 (2) 

where, h(x) = c s x 1 ^ 2 F(x)~ 1 ^ 2 , is half-thickness of the disk and c s is sound speed. 
Now following Mukhopadhyay^, I compute entropy, fi c = M c [M is accretion rate) , 
and energy, E c , at sonic points (x c ) from which are useful to check stability. I 
also compute disk fluid properties solving the set of equations in l|2"jl. are depicted. 



Figure 1. (a) a = (solid curve), a = 0.1 (dotted curve) and a = 0.5 (dashed curve), (b) a = 
(solid curve) and a = 0.998 (dotted curve), (c) a = (solid curve), a = —0.1 (dotted curve) and 
a = —0.5 (dashed curve), (d) a = (solid curve) and a = —0.998 (dotted curve). O and I are 
outer and inner sonic point branches for a = 0. A = 3.3, 7 = 4/3. 

First, following Chakrabartl^ E c — M c analysis has been made in Fig. ^ From 
figure it comes out that with the increase of a, inner sonic points form at more inside 
of inner edge. But the disk looses more entropy at more inside edge, therefore for 
a high a, disk becomes unstable compared to a situation with low or zero a. As 
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the shock formation is directly related to the inner sonic points and its stability, 
for the high 'a' shock becomes unstable or removed (see Chakrabart0 for detail 
about shock in accretion disks). Also from Fig. ^ it is clear that with increase of 
o, transition probability of matter from an outer sonic branch to inner one, with 
increase of entropy, decreases. With increase of counter-rotation, though inner sonic 
points shift to stable high entropy region, the transition probability of matter from 
an outer sonic branch to inner one again decreases and sometimes impossible (Fig. 
^1) . Therefore for both the co and counter -rotating cases the disk and shock tend 
to become unstable. One should keep in mind that nothing is static in Universe 
and a rapidly rotating black hole is most natural. See MukhopadhyajED for detail. 

Figure 2 comes with fluid dynamical results. It directly shows that for moder- 
ately rotating black holes, shock forms at more unstable inside region of disks with 
respect to that of slowly rotating black holes. With more increase of a, any shock 
location approaches to black hole horizon. Therefore for rapidly rotating black holes 
shock disappears completely. Mukhopadhyajd discussed all these in detail. 

Figure 2. (a) Mach no., (b) density and 
(c) temperature in unit of 10 9 for (i) a = 0.1, 
x = 94, Xj = 5.4, A = 3.3 (solid curve); 
(ii) a = 0.5, x = 43, x l = 4.2, A = 2.8 (dotted 
curve) and (iii) a — —0.1, x = 70, x% = 6.324, 
A = 3.3 (dashed curve). Shock invariant quan- 
tity, C, for (d) parameter (i) when x s — 17.18, 
9.2, (e) parameter (ii) when x s = 9.51, 8.35 and 
(f) parameter (iii) when x s = 17.25, 13.73. 
Solid and dotted curves are for inner and outer 
sonic branch. M = 10M Q , M = 1 Eddington 
rate, 7 = 4/3. 



Finally I can say, black hole angular momentum always should be considered, 
before making any inference about global disk properties. It is more important as 
most of black holes are thought to be rotating. Therefore due to the rotation of a 
black hole not only quantitatively but also qualitatively accretion disk properties 
change. 
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